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00 ■ Abstract 

o 

' Previous work on network coding capacity for random wired and wireless networks have focused on the case 

. where the capacities of links in the network are independent. In this paper, we consider a more realistic model, 

5_, ' where wireless networks are modelled by random geometric graphs with interference and noise. In this model, 

Oh. the capacities of links are not independent. By employing coupling and martingale methods, we show that, under 

' mild conditions, the network coding capacity for random wireless networks still exhibits a concentration behavior 

, around the mean value of the minimum cut. 

(N 

I. Introduction 

Traditionally, information flow in networked systems was treated like fluid through pipes, and indepen- 
dent information flows were processed separately. Under this assumption, for a unicast transmission (one 

^ ' source node transfers information to one destination node), the maximum transmission rate is bounded by 
> ■ 

the size of the minimum cut between the source and the destination. This result is known as the Min-Cut 

00 . 

CN ■ Max-Flow Theorem, which was proved by Menger [1], Ford and Fulkerson [2] and Elias et al. [3]. However, 
for a multicast transmission (one source node transfers information to multiple destination nodes), this 
maximum flow rate cannot always be achieved by traditional store-and-forward routing algorithms, even 
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00 

o 



X 



if each source-destination pair has the minimum cut with the same size. That is because in a multicast 
transmission, some links in the network may be shared by the routing paths for different source-destination 
pairs. 

In their seminal paper [4], Ahlswede et al. proposed a network coding scheme, and showed that if 
we allow intermediate nodes to encode their received messages and forward the coded messages to their 
next-hop neighbors, the maximum flow rate can be achieved for mutilcast transmissions. In addition to the 
information theoretic treatment of [4], network coding has also been studied in an algebraic framework 
developed by Koetter and Medard in [5], and a combinatorial framework proposed by Fragouli and 
Soljanin in [6]. Code design for network coding schemes has also attracted intense interest. In [7], Li et 
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al. showed that linear codes are sufficient to achieve the maximum flow rate for a one-source multicast 
transmission. Koetter and Medard, and Jaggi et al. constructed linear multicast codes for network coding 
schemes in [5] and in [8], respectively. The approach of constructing linear codes in a randomized way 
for multicast transmissions was proposed by Ho et al. in [9]. For a detailed review of network coding and 
its applications in many fields, e.g., wireless communication, content distribution, security, please see the 
book by Fragouli and Soljanin [10]. 

In most studies on network coding, network topologies are assumed to be known. In [11], [12], 
Ramamoorthy et al. studied network coding capacity for weighted random graphs and random geometric 
graphs. In the random graph model, each pair of nodes are connected by a bidirectional link with 
probability p < 1 independently [13], [14]. The capacity of each link is assumed to be i.i.d. according 
to some probability distribution. In the random geometric graph model, two nodes are connected to each 
other by a bidirectional link only when their distance is less than or equal to a predefined positive value 
r, the characteristic radius [15]. Each link has a unit capacity. For these two types of random networks, 
the authors showed that the network coding capacity is concentrated at the (weighted) mean degree of 
the graph, i.e., the (weighted) mean number of neighbors of each node. Essentially, the results reveal a 
concentration behavior of the size of the minimum cut between two nodes in random graphs or random 
geometric graphs. Related problems have been studied in the literature, e.g., [16] and references therein. 
In [17], the authors studied a generalized random geometric graph model, where two nodes are connected 
by a bidirectional link with probability 1 if their distance d is less than or equal to tq > 0, and with 
probability p < 1 if tq < d < ri. They obtained similar concentration results. 

The geometric models in [11], [12], [17] assume that a link exists (possibly with a probability) between 
two nodes when the nodes are within each other's transmission range. Although each link has a direction, 
as all links are bidirectional (i.e., the link implies the existence of the link the model in fact 

leads to an undirected graph and considerably simplifies the resulting analysis. In addition, interference 
among wireless terminals was not considered in [1 1], [12], [17]. Nevertheless, in wireless networks, due to 
noise, interference, and the heterogeneity of transmission powers, significantly more sophisticated models 
for link connectivity are needed. For instance, a widely-used model for wireless communication channels 
is the Signal-to-Interference-and-Noise-Ratio (SINR) model [18], [19]. In this paper, we study the capacity, 
i.e., the size of the minimum cut, of random wireless networks under the SINR model. 

Given that network coding capacity with noisy links is in general still an open problem, we assume that 
as long as the SINR jSij of a link (i, j) is greater than or equal to a predefined threshold (3, then node i 
can transmit data at rate R packets/sec to node j without error. That is, links are noise-free once the SINR 
condition is met. In other words, we view network coding as operating on a higher layer in the network 
communication stack, and assume there is an error correcting code at the lower layer which corrects 
errors on the links once the SINR threshold is met. Given this model, each link is indeed directional, and 
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the capacities of different links are not independent. Nevertheless, we will show that under some mild 
conditions, the capacity still has a sharp concentration when the scale of the network is large enough. 

It is worthy mentioning that the capacity we investigate in this paper is different from the one studied 
in [20]-[25]. The latter is referred to as throughput capacity, or transport capacity, for random wireless 
networks with many-to-many transmissions. In other words, it is the maximum achievable averaged rate 
at which each node in the network can transmit (simultaneously with other nodes specified by scheduling 
schemes) to a randomly selected destination node. In contrast, the network coding capacity that we 
study in this paper (as in [11], [12], [17]) is the maximum rate that one source can achieve in a multicast 
transmission, which is determined by the size of the minimum cut between the source and the destinations. 

This paper is organized as follows. In Section II, we describe the random wireless network model. In 
Section III, we study the network coding capacity for a single source and multiple destinations. Specifically, 
we investigate two cases. In the first case, all nodes have the same transmission power, and in the second 
case, the transmission powers are heterogeneous. We use different techniques for these two cases and 
show that the network coding capacity has a concentration behavior in both cases. In Section IV, we 
present relevant simulation results, and finally, we conclude in Section V. 

II. Random Wireless Network Model 
We use the following model for random wireless networks. Assume 

(i) X = {Xi, X2, X„} is a set of n nodes which are independently and uniformly distributed at 
random on the two-dimensional unit torus, where Xj denotes the random location of node i, and n 
is the total number of nodes. 

(ii) Each node i has a transmission power Pi, which follows a probability distribution fp{p), p E 

], where < pmin < 

Here, the existence of a link from node i to node j depends on j's ability to decode the transmitted signal 
from i, which is in turn determined by the Signal-to-Interference-and-Noise-Ratio (SINK) given by 

3.. = (1) 

where Pi is the transmission power of node i, dij is the distance between nodes i and j, and Nq is the 
power of background noise. The parameter 7 is the inverse of the system processing gain. It is equal to 1 
in a narrow-band system and smaller than 1 in a broadband (e.g., CDMA) system. The signal attenuation 
function L(-) is a function of the distance dij = ||Xj — Xj||, where || ■ || is the Euclidean norm, and is 
usually given by L{dij) = cd~°' for some constants c and 2 < a < 4. 

Under the SINR model, the transmitted signal of node i can be decoded at j if and only if [3ij > [3, 
where [3 is some threshold for decoding. In this case, a link (z,j) from i to j is said to exist. Note that 
even if jSij > (3, jSji > (3 may not hold and thus the link (j, i) may not exist. Thus, the graph resulting 
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from the SINR model is in general directed. It is clear that link (z,j) is bidirectional if and only if 
min{/3jj, > (3. Denote by G{X,V,'y) the ensemble of random wireless networks induced by the 
above physical model, where V = {Pi, P2, -P„} represents the set of transmission powers. 
For transmission power P and signal attenuation function L(-), we assume 

(i) Pmin > (3No; 

(ii) Pr(P = prmn) > 0, Pr(P = Pm.ax) > 0, 

(iii) L(x) is continuous and strictly decreasing in x 

for technical and practical reasons. In the remainder of this paper, under different circumstances, we may 
place further constraints on Pj. 

The sum Y^kj^j ^i'^kj) = Y^k^j ^iW^k — Xj||) is a random variable depending on the locations of all 
nodes in the network. Define, for all j = 1, ...,n, 

J(j)=5^L(4,), (2) 

k^j 

J(j) = 5^Pfci^(4,). (3) 

ki^j 

To study the asymptotic network capacity, we will let the number of nodes n go to infinity. Since the 
region is fixed, this corresponds to a dense network model [15], [20]. Another widely used model is the 
extended network model [26], [27], in which the number of nodes and the area of the region both go 
to infinity while the ratio between them — the density of the network, is kept constant. Both models are 
widely used in the literature. We will focus on the former one in this paper. 

III. Network Coding Capacity for Single Source Transmission 
A. Capacity of a Cut 

Let Cij be the capacity of a link We will specify the form of Cij later for different scenarios. 

Consider a single-source multiple-destination transmission problem. Let s be the source node. Suppose 
there are / destination nodes, ti, ...ti, and m relay nodes, Ui, Denote the set of the destination nodes 
and relay nodes by T and TZ, respectively. Note that {s},T and TZ are all subsets of X. In this paper, 
we always assume that there are no direct links between the source and its destinations. Fig. 1 gives an 
example of single-source single-destination transmission. 

Let the capacity of the link from the source s to each relay node Ui be Cgi, i = 1, ...,m, the capacity 
from relay node Uj to another relay node uj be Cij, i ^ j,i = 1, m,j = 1, m, and the capacity from 
each relay node Ui to each destination node tj be Cu-, i = 1, m, j = 1, /. Unlike random geometric 
graph models studied in [11], [12], [17], the capacities in our model are not symmetric (i.e., Cij ^ Cji) 
nor independent in general. 
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Fig. 1. Single-source single-destination transmission in directed SINR graplis 

In our SINR wireless network model, there are two sources of randomness: one is the random location 
of each node and the other is the random transmission power of each node. We use Ex and Ep to denote 
the expectation operation with respect to each probability measure, respectively. 

Let C be the expected capacity of a link (i, j), defined as 

C = ExEp[Cij\ 

Q.dFpJr), (4) 



where F|3^.{■) is the c.d.f. of f3ij, which is determined by fp{-), the distribution of X, and the path-loss 
function L{-). 

Now define an s-t-cut of size k for a given source s and destination t G T as a partition of the relay 
nodes into two subsets and Vj^, such that \Vk\ = k, \Vi^\ = m — k, VkV^V^ = TZ and fl V{: = 0. An 
example of an s-t-cut is shown in Fig. 2. Let 

then Ck is the capacity of the corresponding s-t-cut. Although Ck is a sum of dependent but identically 
distributed random variables, we still have 

E[Ck] = ExEp[Ck] 

= ^ ExEp[Cs^] ExEp[Cj,] + J2 ExEp[Cjt] 

= [m + k{m - k)]C, (6) 

and consequently E[Ck] = E[Cn,^k] for A; = 0, 1, ...,m, and E[Co] < E[Ci] < ■ ■ ■ < E[Cim/2-]]- 

To show that the capacity of any source-destination pair concentrates at some value, we will first show 
that for such a source-destination pair, the capacity of any s-t-cut of size k concentrates at its mean value. 
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Fig. 2. An s-t-CMt for the single-source single-destination transmission in directed SINR graphs 



Similar results were proved in [11], [12], [17], where the capacities of the links that originate from the 
same node are i.i.d. Nevertheless, the methods used in [11], [12], [17] do not apply here, since in the 
SINR model, Ck is a sum of dependent random link capacities. Instead, we employ coupling, martingale 
methods and Azuma's inequality [15], [28] to solve the problem. 

Note that when 7 = 0, i.e., there is no interference in the network, the capacities Csi for i = k + 1, m, 
are mutually independent; so are the capacities Cij for any fixed i = 1, ...,k with j = k + 1, ...,m or t. 
In this case, although the link capacities are still asymmetric, Xlwiev^? ^si and Yju,ev^u{t} ^i* j ^ 
become sums of independent random variables. Thus we can apply methods similar to those used in [11], 
[12], [17] to obtain the same concentration results. 

B. Constant Transmission Power 

Consider the scenario when all nodes transmit with a constant power Pq and denote the corresponding 
model by G{X, Pq,^). In this case, the SINR of link (i, j) can be rewritten as 



A- 



L{dij) 



(7) 



No/Po + lJ{j)-im,y 
Assume that when Pij > (3, the link (z,j) has capacity R, i.e., node i can transmit data at rate R 
packets/sec to node j without error. Then, we can define as 

~ \ A, < (5. 

Note that when the wireless channel is an additive Gaussian channel, the capacity of link (i, j) is [29] 

Ai < 
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Our results in this subsection do not depend on any particular form of Cy when [3ij > (3. Nevertheless, 
since we consider the application of network coding, it would be more appropriate to focus on dS]), rather 
than dll). 

Note that jSij and thus Cij, are determined by L{dij) and J{j). Because of the i.i.d. distribution of the 
Xj's, given Xj, the dij's are independent for all i ^ j. Given node j, let 

E[L] ^ E,,Md^J)], (10) 

then 



{n-l)E[L]^E[J]. (11) 



E[jij)] = E 

Since our model is a dense network model and the area of the region is fixed, E\L\ = E[L(dij)] 
is a constant and E[J] = [n — \)E[L\ scales with n. For different j's, it is clear that J{jys are not 
independent. However, they have the same sharp concentration behavior in large-scale wireless networks. 
This is established in the following lemma. 

Lemma 1: Suppose there are n nodes in the network, then 

Ft{J{j) < {1 - e,)E[J]) = O ( ^) , (12) 

and 



Fr{J{j)>{l + e[)E[J])=0[\], (13) 



41nn ntiA / 61nn 



for all J = 1, 2, n, where ei = and - ^ (^^n^ ■ 

Proof: Given any node j, because J(j) = ^i'^ij)^ and L{dij) are i.i.d. for all i ^ j, by the 

Chernoff bound [28], [30], we have 



PT{JU)<{l-e,)E[J])<exp 



m4 

2 



. {n-l)E[L]el. 



and 

Pr(J(j) > {l + e[)E[J])<exp 



me? 

3 



, (n - l)E\L]e\^ , 
exp<;-^ ^"^^i^- ^^^^ 



Substituting ei = ^J (^.^1^^e[l] = y {n-i)E[L] ^ '^^ obtain ^ and ([T3l) . 

respectively. □ 
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Lemma [T] shows that when the network is large, i.e., n is sufficiently large, J(j) concentrates at 
7(n — 1)E[L] = Q{n)E[L]. The reason for this is the uniform distribution of the nodes. 

Now define two other types of SINR models G'{X , Pq,'j) and G"{X , Pq,^) which are coupled with 
G{X, Pq, 7) such that they have the same point process X and constant power Pq. Let the SINR of link 

in G'(A',Po,7) and G"{X,Po,i) be 



~ iVo/Po + (l + e'i)7^[J] - lL{d,^) 



A', = , , \ 'Lr. , ... ^ (16) 

and 



B". = (17) 



respectively. 

Let Glj and C-^ be the capacity of link in G'{X, Pq,^) and G"{X, Pq,^), respectively. Since 

ei ^ and e'^ ^ as n ^ 00, G-j and G-'j are asymptotically equal to Gij. 

The following lemma establishes a concentration result for Gk with constant transmission power by 
coupling methods. 

Lemma 2: For any < e < 1, the capacity of an s-t-cut of size k,k = 0, 1, m, satisfies 

Ft{G, < (1 - e)E[G',]) < exp ( -^^^] ( 1 - O ( -)) + O ( -) (18) 



where -E'[C[,] = [m + k{m — k)\C' and C' is the average link capacity in G\X, Pq, 7). Moreover, 

Pr(C. > (1 + e)E[G'^]) < exp [-^^^] (l - O + O (19) 

where -E[C[f] = [m + A;(m — A;)]^" and C" is the average link capacity in G"{X, Pq, 7). 

Proof: By Lemma [B for all j, {J{j) > (1 - ei)^[J]} and {J{j) < (1 + e;)E[J]} are both increasing 
eventslllBy the FKG inequality [15], [26], [30], we have 

(n \ n 

f]{J{j)>il-e^)E[J]}] > l[^r{J{j)>{l-er)E[J]) 



1-0 



1 



1-0 1-], (20) 



n 



'in the context of graph theory, an event A is called increasing if Ia(G) < Ia{G') whenever graph G is a subgraph of G', where Ia is 
the indicator function of A. An event A is called decreasing if A'^ is increasing. For details, please see [15], [26], [30]. 
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where the first equality is due to Lemma [IJ Similarly, 

/ n \ n 



\j=i / j=i 

' 1 



1-0 



'n? 



l-O(-). (21) 



n 



Inequalities (|20l) and (|2T]) imply that 

Pr(Cfc <Cl)>l-0 , (22) 

Pr(Cfc >C'^)>l-o(-\. (23) 



Since 



and 



Pr(C, < (1 - e)E[C',]) < 1 - Pr(C, > Pr(Cl, > (1 - e)E[C',]), 



Pr(C, > (1 + e)i?K]) < 1 - Pr(C, < C^') Pr(C[' < (1 + e)E[C'^]), 
in order to show (fTSl ) and ( fT9l ). it suffices to show 

Pr(C^ < (1 - e)E[C',]) < exp {-^^^fi^} , (24) 

and 

Pr(C^' > (1 + 6)EK1) < exp {-^^ifl^} • (25) 

In G'{X, Po, 7) and G"{X, Pq, 7), the SINK of Unk (i, j) is given by (dH) and respectively. Because 
the dij's for a given i are independent, by applying the Chemoff bounds, we obtain (|24|) and (l25l) . □ 

Since C^'^ and C-^ are asymptotically equal to C.^, -E'[C[] and -EfC^f] are asymptotically equal to E[Ck]. 
Consequently, Lemma |2] shows that Ck concentrates at E[Ck] asymptotically almost surelyH 
Now, let Cs^t be the minimum cut capacity among all s-t-cuts, i.e., 

Cs,t = min Ck- (26) 

0<k<m 

For one source and multiple destinations, the capacity of network coding depends on the minimum cut 
between the source and the destinations [11], [12], [17]. Therefore, for the given source node s and the 
sets of destination nodes T = {ti,...,ti} and relay nodes TZ = {ui, ...,Um}, define the network coding 
capacity as 

Cs,r = minC^f. (27) 

tGT 

^An event is said to be asymptotic almost sure (abbreviated a.a.s.) if it occurs with a probability converging to 1 as n — > oo. 
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In the following, we show that when the number of relay nodes m is sufficiently large, the network 
coding capacity Cs^t concentrates at -^[Co] = mC with high probability. 

Theorem 3: When n is sufficiently large, with high probability, the network coding capacity r 
satisfies 

' (28) 



Pr(C,r>(l-6M^^o]) = l-o(^) 



where = -^^f^ for a > and E[Co] = mC. 

Proof: Since the Ci/s are asymptotically equal to the C-^'s, in order to show (|28]) , it is equivalent to 
show 

Pr(C.,r > (1 - e'MC'o]) = 1-0 . 
Since E[Cl] > for any k = 1, ...,m, 

Pr(C,, < (1 - e'jE[C'o]) < FT{Cs,t < (1 - e'jE[C'„]), 

for any t E T, where k' is the size of the minimum s-t-cut. By (fTSi) of Lemma [21 we have for sufficiently 
large n, 

e'^[m + k'{m-k')]C'' 



Pt{C,, < [1 - e'jE[Ci,]) < exp 

< exp 



2 



2 



By choosing e'^ = y ^^[^'p since C and C are asymptotically equal, we have for any t E T, 

Pr(a, < (1 - e'^)E[C',]) = O . 

By the union bound, we have 

PT{C,,r<{l-e'^)E[C'o]) < Y^PriC,, < {1 - e'jE[C',]) 

I 



o 



□ 



Theorem 4: When n is sufficiently large, with high probability, the network coding capacity Cg^r 
satisfies 

Pr(C.,r < (1 + el)E[C,]) = 1-0 , (29) 



where e'^ = for a > and E[Co] = mC. 
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Proof: Since the C^j's are asymptotically equal to the C^^'s, in order to show (|29l ). it is equivalent to 
show 



Pr(C.,r<(l + 6'^)i?ra) = l-0(^-L). 



To show this, it is sufficient to consider a particular cut for a source-destination pair, e.g., an s-t-cut with 
capacity YlT=i ^si separating the source s from all the other nodes. 

Pr(C.,^ > (1 + e:)i?[Q) < Pr (C,, > (1 + 6'^)i?[Q) 

(m 

= Pr(Co>(l + e:)E[Q) 



< exp 

- of- 

where the last inequality follows from (fT9l) . □ 



C. Heterogeneous Transmission Powers 

In this subsection, we consider the case where the transmission power of each node is randomly chosen 
rather than being constant. We continue to assume that the capacity of a link is a constant R, which 
is independent of the SINR Pij, when Pij > (3. In this case, Pij can be rewritten as 



PiL{dij ^ 



N0 + lEk^^,Pkm,) 

(30) 



No + lI{J)-lPMd^J)' 

Because the -Pj's and Xj's are both i.i.d., using the same technique as that for Lemma [H we can prove 
the following lemma. 

Lemma 5: Suppose there are n nodes in the network, then 

Pr{I{j)<{l-e2)E[I])=o(^], (31) 

and 



Pr(J(j)>(l + 4)^[/])=0 - , (32) 



41nn ^' / 61nn 



for all J = 1, 2, n, where 63 = ^J (n^i)E[P]E[L] - y ^n^i)E[P]E[LV 

Even though we have concentration results for /(j), we cannot employ the same coupling methods 
as in the previous section. In G'{X,Po,^) (or G" {X , Pq, -y)) as described in Section III-B, the C^'^'s (or 
C"j's) are independent for all j ^ i for a given i. In our new case, however, this independence does not 
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hold because all the Cj/s depend on the transmission power Pi. To deal with this dependence, we use 
martingale methods and Azuma's inequality to solve our problem. 

Theorem 6 (Azuma's Inequality [30]): Let Zq, Zi, he a martingale sequence such that for each i = 
1,2,...,, 

\Zi — Zi_i\ < Cj 

almost surely, where Cj may depend on i. Then for all n > and any A > 0, 



2ELiC? 



Pr(Z„ >Zo + X)< exp <| , \ , (33) 



and 

A2 



Pr(Z„ < Zo - A) < exp . (34) 

Proof: Please see e.g. [30]. □ 

To use Azuma's inequality, we need to construct a martingale. A common approach to obtain a martin- 
gale from a sequence of (not necessarily independent) random variables is to construct a Doob sequence. 
Specifically, suppose we have a sequence of random variables Yi,Y2, ...,Yn, which are not necessarily 
independent. Let S = J2^=i ^ define a new sequence of random variables {Zi : i = 0,1, ...,n} by: 

Zq = E[S] 

Z^ = Ey,^„...,y„[S\Yu...,YI i = l,2,...,n. ^^^^ 
Then {Zt : i = 0, 1, n} is a martingale and Z„ = 5*. 

If we are able to upper bound the difference \Zi — for all i by some constant, then we can apply 
Azuma's inequality to obtain a bound on the tail probability. For example, if the Fj's are independent, a 
simple upper bound for \Zi — Zi^i\ is any upper bound on \Yi\. However, if the Kj's are dependent, we 
need to further understand the properties of the Fj's to bound \Zi — We approach our problem by 

following this idea and using the next corollary. 

Lemma 7: For n > 1, given a sequence of random variables Yi,Y2, ■■■,Yn, which are not necessarily 
independent, let S = Yi. If for any yi, y'- E Di, where Di is the support of Yi, 

\E[S\Y,, Y, = y,] - E[S\Y,, Y, = y'^] \ < q, 

almost surely, where Cj may depend on i, then for any A > 0, 



Pt{S > E[S] + A) < exp { , \ , (36) 



and 

Ft{S < E[S] - A) < exp {-^^^ } ■ (37) 
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Proof: We prove this corollary for the case of discrete random variables. For continuous random 
variables, the proof is similar. 

Define a Doob sequence with respect to {Yi : z = 1, as in (l35l) . To simplify notation, we will 

write -E'y,+i,...,y„[5'|Yi, as E[S\Yi, ...,Yi] when there is no ambiguity. 

By the total conditional probability theorem, we have 



and 



Z, = E[S\Yi,...,Y^ 

= J2E[S\Y^,...,Yi\PT{Y, = y\Y,,...,Y, 



i-i) 



Therefore, 

\Zi — Zi^i\ 



< Yl \E[S\Y,, F.] - E[S\Y,, = y]| Pr(F, = 



y6A 
Ci. 



Since {Zi : i = 0, is a martingale with bounded difference of \Zi — we can apply 

Azuma's inequality to obtain (l36l) and (1371) . □ 



Now consider G'{X,V,^) and G"{X,V,'y) coupled with G{X,V,'y) such that they have the same 
point process X and powers V, where the SINR's of link (z, j) in G'{X,V,'y) and G"{X,V,'y) are 

PiL{d 

iVo + (1 + 4)t^/] - 7^.^(c?«j) 

and 



A'. = . , . (38) 



nil ^ ^ ^ ''^^ (39) 

No + il-e,hE[I]-^P,m,y ^ ^ 

respectively. 

Let G'ij and be the capacity of link (ij) in G"(A',P,7) and G"{X,V,-f), respectively. Then, 
and G'/j are asymptotically equal to Gij. 

Assume that there exist positive solutions r'„j„ > 0, r^^^^ > 0, r^j„ > and r^„^, > for the equations 



No + ^{l + e',)E[I]-^P:^Mr'^J 
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and 



respectively. That is 



No + 7(1 - 








iVo + 7(l 


-e2)E[I] - -fp^inL{r'^-J 







iVo + 7(l-e2)i?[/]-7Pn^a.i^(C 



r' . 

mm 



P No + ^il + e',)E[I] 



1 + 7/5 


Pmin 








)E[I] 


1 + 7/5 


Pmax 






No + 7(1 - 




1 + 7/? 


Pmin 






No + 7(1 - £2, 





Pmax 



,1 + 7/5 

Since L(-) is continuous and strictly decreasing, r^j„, r'^ax^ ^min ^^^d r^^^ are all unique. In G'{X, V, 7) 
(G"{X,V,^)), any node inside the circle centered at with radius r'^j„ (r'min) connected to node i 
by a bidirectional link, while any node outside the circle centered at Xj with radius r^^^ i^max^ '^ot 
connected to node i. 

Let ^(Xj, r^j„, r'^ax) ^(Xj, r^j„, ^^ma^,) be the two annuli with inner radius r^j„ and outer radius 



and inner radius r;;^^ and outer radius C^^,, respectively. Denote by A^(Ci„, C^^.) and A^^^^^, r 



max J 



the number of nodes in ^(Xj, 



and^(X„r;;^ 



), respectively. It is clear that iV(r;^.„, r^^^^) 



and N{r'^^ 



mm^ max. 



are binomially distributed with mean mi{r 



/2 



max mini "^^^^ V max mm. 



JI2 



, respectively. 



Now suppose the signal attenuation function L{x) 
Then, 



cx 



c(l + 7/5)p, 



for some constants c > and 2 < a < 4. 



1/0 



and 



„'2 



^max ^min 



/?[iVo + 7(l + 4)^W] 
c(l + 7/5)p 

max 

/?[iVo + 7(l + 4)^W] 
P[No + iil-e2)E[I]] 

C(l + -i(3)pmax 

f3[No + l{l-e2)E[I]] 

E {Pmin 1 Pmax 



1/a 



1/a 



1/a 



^"2 _ ^"2 

max min 



[iVo + 7(l + 4)^[/]]'/° 

^ iPmin ; Pmax ) 

[iVo + 7(l-e2)^[/]]2/'^ 



(40) 
(41) 
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where 

-D\J>min, Pmax ) — [Pmax Pmin. 



P 

From (|40|) and ((4T|) . we can see that both mr{r'j^^^ — r'^in) ^^^d n7r(r^^^^ — r^^j„) scale with n as — ^ 



(42) 



max minJ "i^" V maa; minJ '^^^^^ vviiii lu cio 5 

since E[I] = [n — 1)E[P]E[I]S scales linearly with n. Now assume that there exists a constant 77 > 
independent of n such that 

N{r'^in,r'^ax) and A^lr;;^^, C^J < ^ (43) 

a.a.s. This assumption effectively puts a constraint on the transmission power. For example, if we choose 
?7 = 1, then the transmission power P must scale with n so that Pmax = 0(n^~°/^). Note that r^-„ and 
'^'max arc asymptotically equal to r^-„ and r'^^x, rcspectively. 

The following lemma establishes a concentration result for Ck with heterogeneous transmission power 
and constant capacity by coupling methods and Azuma's inequality. 



Lemma 8: For any < e < 1, when n is sufficiently large and (|43l) is guaranteed, with high probability, 
the capacity of an s-t-cut of size A;, /c = 0, 1, m, satisfies 

Pr(a<(l-.Wg])<e.p{ J'"^'''-f'l^°'^ }, (44) 

where E\C'f^ = [m + k{m — k)]C' and C" is the average link capacity in G'{X,V,'j). Moreover, 

where E[C'^] = [m + k{m - k)]C" and C" is the average link capacity in G"{X, V, 7). 

Proof: By Lemma |5l for all j, (1 - e2)E[I] < < (1 + e'^)E[I] holds a.a.s. As in the proof for 
Lemma [2l we have that is lower bounded by C'f., and upper bounded by C'^, with probability at least 
1 — 0{l/n). Hence, in order to show (f44)) and (l45l) . it suffices to show 

p.(c; < (1 - .)E|cy) < exp |_ l- + M.^-^]c?'V | ^^^^ 

and 

P..(q' > (1 + < exp{- l" + ^<-f)l^"-' } , (47) 

To show (061), we use martingale methods. Let Yi = C^(fc_^^), F2 = C'^(k+2)^ ^^-fc = '^L' ^m-fc+i = 

^l(fe+l)' ^^-A:+2 = C*i(fc_|_2)! ■^-fc+fc(m-fc) = and Ym~k+k{m-k)+l = C*lt' ^-A:+fc{m-fc)+2 = 

C'2i., ■..,Ym-k+k{m-k)+k = C*L- Define a Doob sequence {Zi : z = 0, ...,m + A;(m — k)} with respect 
to {Fj, i = 1, 2, m + fc(m — fc)} as 

Zo = E[C',] 

Z, = E[C'k\Y,,...,Yi\, i = 1,2,.. .,m + k{m-k). 
Then {Zi : i = 0, ...,m + k{m — /c)} is a martingale and Zm+k{m.~k) = C'^. 
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Since when i ^ u and j ^ v, C-^ is independent of C'^^„ dependence exists only among C(j's, j ^ i, 
for a given i. However, the distances dij are independent for all j ^ i for given i. Hence the difference 
between E[C'f^\Yi, ...,Yi_i,Yi = yi] and £'[C^|Yi, = yl] depends only on those F^'s which are 
incident on the same source node as ¥[. When rf^ < r,^j„, C'^j = R, and when dij > r'^^ax^ C[,j = 0. 
Moreover, the number of nodes within the annulus A(X.i, r'^i^, r'^ax) upper bounded by the constant i] 
a.a.s. Therefore, 

\E[C',\Y,,...,Yl.,,Yl = yl]-E[C',\Y^,...,Yl_,,Yl = yl]\ 

< vR 

a.a.s., where yi and yl are either or R. 

Applying the result of Lemma Ul we have (|46l) . In the same manner, we can show that (|47] ) holds. □ 



In the following, we show that when the number of relay nodes m is sufficiently large, the network 
coding capacity r concentrates at E[Co] = mC with high probability. The proofs are based on Lemma 
[8] and very similar to those for Theorem [3] and Theorem |4l and provided here for completeness. 

Theorem 9: When n is sufficiently large, with high probability, the network coding capacity Cg^r 
satisfies 

Pr(C7,,T > (1 - e^)E[Co]) = 1-0 ( ^) , (48) 

where = ^^^^^ V2am Inm for a > and £"[6*0] = mC. 

Proof: Since the C^/s are asymptotically equal to the C^^'s, in order to show (l48l) . it is equivalent to 
show 

Pr{Cs,r > (1 - e^)E[C',]) = 1-0 (^-^^ . 
Since E[C'f^] > E[C^] for any A; = 1, ...,m, 

Pr(C,,< < (1 - e^)E[C',]) < Pr(C,,t < (1 - e^)E[C',,]), 

for any t E T, where k' is the size of the minimum s-t-cut. By (|44)) of Lemma [8l we have 

€l[m + k'{m-k')]C' 



PT{Cs,t < {I - ea)E[Ci,]) < exp 



< 



2(r/ + l)2i?2 



exp 



2(r/ + l)2i?2 



By choosing = ^ ^Q^"^ ^ a > 0, since C and C are asymptotically equal, for any t G T, 



Pr(C., < (1 - 6jE[C^]) = O . 
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By the union bound, we have 

Vi{Cs,r<{l-eo)E[C'^]) < ^ Pr(C,,< < (1 - e,)i?[C^]) 




□ 



Theorem 10: When n is sufficiently large, with high probability, the network coding capacity Cs,t 
satisfies 

Pr(C,,r < (1 + e^)E[Co]) = 1-0 ( ^] , (49) 



where = , \j2am Inm for a > and -EfCo] = vnC . 



Proof: Since the Cjj's are asymptotically equal to C'-^'s, in order to show (|49l) , it is equivalent to show 



Pr(C,,^ < (1 + e^)E\Cl\) = 1-0 . 



To show this, it is sufficient to consider a particular cut for a pair of the source and one destination, for 
instance, an s-t-cut separating the source s from all the other nodes. 

Pr(C,,T>(l + e,)E[C^]) < Pr(C,,t>(l + e„)E[Q) 

/ m 

< PrK^C,,>(l + e„)i?[Q 



, i=l 



= Pr(Co>(l + e,)Era) 
< exp 

= O 



'2(r/ + l)2i?2 

1 

where the last inequality follows from (l45l) of Lemma [H □ 



IV. Simulation Studies 

In this section, we present some simulation results on the SINR model and the network coding capacity. 
Fig. 3 and Fig. 4 show simulation results on interference and cut capacity in G{X, Pq, 7), where n = 2000, 
L{x) = ^^^x~^, Nq = 0.02, P = 0.2 and 7 = 0.02, and every node transmits with constant power 
Pq = 0.01. Fig. 5 and Fig. 6 show simulation results on the interference and the cut capacity in G{X, V, 7), 
where n = 2000, L{x) = A'^o = 0.02, (3 = 0.2 and 7 = 0.02, and every node transmits with 

power P uniformly randomly distributed over [0.01, 0.02]. The results confirm the concentration behavior 
of the interference and the cut capacity. 
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Interlerence Level Capacity ot s-t-Cut ot Size 50 (1<=50) 



(a) Interference J{j) (b) Cut capacity Ck 

Fig. 3. Interference at eacii node, and capacity of random s-t-cut of size fc = 50 in G{X , Po,-y) 



% 50 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 500 600 700 800 900 1000 1100 

Interference Level Capacity ot s-t-Cut ot Size 50 (k=50) 



(a) Interference J{j) (b) Cut capacity Ck 

Fig. 4. Interference at each node, and capacity of random s-t-cut of size fc = 50 in G{X jVyj) 

V. Conclusions 

In this paper, we studied the network coding capacity for random wireless networks under a SINR 
model, where the network is modelled by the graph G{X ,V,'j). Previous work on the network coding 
capacity for random wired/wireless networks are based on the assumption that the capacities of links are 
independent. In the SINR model, however, the capacities of links are not independent due to noise and 
interference. We investigated two scenarios. In the first case G{X, Pq,^), we assumed all nodes transmit 
with a constant power Pq. To study the network coding capacity in G{X, Pq, 7), we coupled G{X, Pq, 7) 
with two other models G'{X, Pq, 7) and G"{X, Pq, 7), which have the same point process X and constant 
power Pq, but different thresholds. We showed that the network coding capacity for G{X, Pq, 7) is upper 
and lower bounded by those for G'{X,Po,'y) and G"{X , Pq^'j). By proving that the network coding 
capacities for G'{X , Po,'j) and G"{X , Pq,^) concentrate on the same value asymptotically, we showed 
that when the size of the network is sufficiently large, the network coding capacity for G{X, Pq, 7) exhibits 
a concentration behavior around the mean value of the minimum cut. In the second case G{X ,V,'j), we 
assumed each node transmits with a random power drawn from some distribution. Since coupling methods 
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could not be applied in this general case, we used martingale techniques to deal with dependence between 
link capacities, and showed that under some mild conditions, the network coding capacity also exhibits a 
concentration behavior. The results obtained are important for understanding network coding performance 
in random wireless networks under the SINR model. In addition, the methods used in this paper provide 
useful techniques for studying properties of random wireless networks under the SINR model. 
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